Statics and dynamics of the Lebwohl-Lasher 
model in the Bethe approximation 



N.S. Skantzos 1 

Instituut voor Theoretische Fysica, Celestijnenlaan 200D, Katholieke Universiteit 

Leuven, B-3001, Belgium 

and 

J.P.L. Hatchett 2 

| Hymans Robertson LLP, One London Wall, London EC2Y 5EA, UK 
Abstract 

We study the Lebwohl-Lasher model for systems in which spin are arranged on 
random graph lattices. At equilibrium our analysis follows the theory of spin-systems 
on random graphs which allows us to derive exact bifurcation conditions for the 
phase diagram. We also study the dynamics of this model using a variant of the 
dynamical replica theory. Our results are tested against simulations. 
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1 Introduction 



Some physical systems are known to exhibit remarkably sharp, and yet continuous, 
transitions between the ordered and paramagnetic phase. The most prominent ex- 
ample of these are liquid crystals, systems which combine order like that found 
in solids with fluidity like that of liquids. Liquid crystals, generally modeled as 
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systems of 'hard-rods', typically exhibit a transition between a phase with no ori- 
entational or translational order and a phase where an ordered structure appears. 
The nature of this transition (continuous vs discontinuous) is an important aspect 
in practical applications. 

One of the most successful models that is able to capture the main character- 
istics of this transition was introduced by Lebwohl and Lasher in 1973 [1]. In 
this model the microscopic degrees of freedom & take real values in the interval 
[0, 27r), representing the orientation of rod i relative to some fixed reference point. 
To model the sharp transition, the coupling energy between any pair of rods is 
taken to have the shape of a deep and narrow well. This energy is taken to be 
€ij{(f>i — (f)j) = —JL p (cos(0j — 4>j)) for any pair of rods where J represents 

the strength of the interaction and the function L p (x) denotes a p-th order Leg- 
endre polynomial. In the context of 'hard-rods', p is taken to be even to enforce 
invariance of the energy under the transformation — > + n. The value of 
p plays a crucial role in the nature of the phase transition [2,3]. Generally, this 
model has been the subject of a significant amount of research (see for instance 
[2,3,4,5,6,7,8,9] and references therein). 

In this paper we study the Lebwohl-Lasher model on a sparse graph motivated 
both by the interesting nature of the transition and the relative scarcity of ana- 
lytic results. Our approach can be viewed as the Bethe approximation to the finite 
dimensional problem, whereby for every site there is an explicit local neighbor- 
hood. At equilibrium our analysis follows the finite connectivity (as opposed to 
fully-connected mean field) theory [10] as developed for real- valued spin systems 
[11,12]. In the thermodynamic limit we are able to solve this model exactly and 
derive expressions describing the bifurcation lines in the phase diagram. Numerical 
evaluation of the order parameters agrees well with the results of the bifurcation 
analysis and shows that the phase transition, although sharp, is second-order. 
Comparison of the above analytic results with simulation (Langevin) expreriments 
shows excellent agreement. 

We have also studied the dynamics of this model. Here, in contrast to the ther- 
modynamic analysis which follows a well-studied framework, the terrain of the 
finite-connectivity dynamics is much less explored. Some recent advances are the 
papers [13,14,15]. In the present study we have chosen to extend the results of [13] 
in order to account for the continuous nature of the spin variables. As in the dy- 
namical replica theory of [16] on which our analysis is based, we arrived at closed 
equations by assuming equipartitioning of the microscopic state probability within 
the observable subshells. The resulting analytic description is in good agreement 
with Langevin simulation for small and very large times, but not sufficiently accu- 
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rate for intermediate times. This is an artifact of the equipartitioning ansatz and 
the truncation of the set of observables to a relatively small (and computationally 
tractable) set. 



2 Model Definitions 

The model consists of N microscopic variables <fi = (fa, . . . , <Pn) representing the 
angular phase of oscillators relative to a fixed frame of reference. At equilibrium 
the system is described by the Hamiltonian 

H(<t>) = -J J2 L p (cos((f)i - (f)j - uj i:j )) (1) 
(i,j)eg N 

where Qn = (Vat, £n) denotes a graph instance of a set of vertices = {1, . . . , N} 
and edges Sn = £ Vat}. The function L p (x) denotes a p-th order Leg- 

endre polynomial. The thermal variables take values from the interval fa G [0, 27r) 
for alH = 1, . . . , N. We also introduced the 'disorder' angle variables Uij G [0, 2n) 
which represent locally preferred orientations between sites % and j. In a physical 
scenario these are induced by the presence of impurities, cavities, or other hetero- 
geneities. We will assume that the distribution of these variables is for all i < j 
given by 

K (Vij) = \ + \ <L« ,-ffl (2) 

for some uo G [0, 2ir). A more general treatment is straightforward provided that 
the detailed balance condition is met. For each spin in the system we assign a local 
neighborhood, i.e. a set of sites to which the spin is connected. Let us for simplicity 
abbreviate this set using di = {£\(£,i) G Sn}, with \di\ representing the size of the 
neighbourhood of i. Our graph is characterized by the 'degree' distribution 

i 

Moments of this distribution will be denoted (k m ) = J2k>oP(k)k m - We draw graphs 
uniformly from the ensemble of all graphs that have this given degree distribution, 
i.e. we use the configuration model [17]. In the thermodynamic limit this leads to 
interactions on a sparse graph for which the Bethe approximation is exact. 
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3 Equilibrium analysis 

3. 1 The cavity formalism on a single instance 

To evaluate thermodynamic properties and the phase diagram of this system we 
have two main analytic tools at hand: the replica [18] and the cavity method [10]. 
Both of these have been used widely in a variety of settings, both in infinite- 
and finite-connectivity systems. Although the underlying philosophy of the two 
methods is different they deliver identical results. In this paper we will choose to 
follow the cavity method as developed for real- valued variables in [11]. 

The starting point of this method is to consider the marginal probability of finding 
a spin at site i in state fa. This follows from 



where Zi is the appropriate normalisation constant. (<!>&) is the joint probability 
of finding the neighbours of i in state 0^ = (0^, . . . , (f>j, gi .) in the absence of site 
%. To understand the origin of this expression we imagine an iterative process that 
creates an exact tree structure in which each spin is connected to k neighbours 
apart from the site at the top of the tree which has k — 1 neighbours. At each 
step of this process we add a new spin to the system by bringing together k — 1 
independent branches and connecting the top site of each of these to the new site. 
Since, the k — 1 branches influence each other only through the new site that has 
been added, we can view the probability P^((j>di) as the a priori probability of 
finding 0^ before the addition of the new spin, whereas the Boltzmann factor 
accounts for the energetic cost of the merging. Since before the addition of the 
new site into the system, the k — 1 branches were independent we can factorize 



This equation is the Bethe approximation and is exact, by construction, on trees. 
For random graph structures where loops, although predominantly long (e.g. scal- 
ing as log N on Erdos-Reyni random graphs), do exist, (5) is approximately correct. 

Equation (4) is not closed as it relates the true probability distribution with the 
cavity probability distribution. To close it we simply remove a neighbour of i which 





(5) 
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leads to 



The two types of distributions that appear above can be dealt with either (i) by 
an explicit parametrization, namely setting P^^(0j) — > P(</>j|{/4^}) (and similarly 
for P(4>i)) where {/4^} play the role of cavity fields [11,12] or (ii) by using a sim- 
ple histogram for each of the distibutions. Note that since spins are continuous 
variables, we are in principle required to introduce an infinite number of cavity 
fields. For all practical purposes however, one truncates the number of fields, while 
in some cases appropriate choices for the parameterization can be made which 
ensure that the impact of the truncation is not significant [12]. The histogram 
method of course requires advanced computational power since one needs to allo- 
cate memory space for a relatively large number of bins, but this approach has the 
advantage of working directly with the distributions. 

Once the stationary values of (6) are known for alH = 1, . . . , N and j G di we 
can evaluate the true probability function (4) and subsequently observables in the 
system, e.g. the 'magnetizations' 

H-i^e/^W 008 ^] (T) 

\m,) "~°° N < 1 \M<P.) j 

and similarly for other order parameters. 



3.2 The cavity formalism in the graph ensemble 



To perform a bifurcation analysis of this system it is convenient to work with a 
graph ensemble, i.e. the set of all instance graphs Gn = {Gn} with the given 
degree distribution (3). 

We now consider the population of the cavity probability distributions in the graph 
for all sites. This defines a functional density in the following way 

i I I j£di 

where by 6(p) we mean a functional delta distribution in the sense that w[f] = 
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/ dgw[g] 6(F)\g(x)—f(x)]. Similarly, for the density of true probability distributions 



W tIue [{P}} = ton -£<f (F) [P(0) -P,(0)] 



(9) 



We can now convert (4) and (6) into self-consistent equations for (8) and (9) for 
the ensemble of graphs with degree distribution p(k) 



k-l 



l\[{dP e }W[{P e }} doj e K(u e )\ 



e=i 



x 5 



(F) 



-i k—1 

PW II [W W) eP JL ^ cos( ^'-^ 



(10) 



where Z is the appropriate normalization constant. The prefactor p(k)k/{k) ex- 
presses the fact that the probability that any given site is chosen to become a 
cavity site is proportional to the number of bonds it has. For the density of true 
probability distributions we have 



fc>0 

x 5, 



W tIUC [{P}] = J2p(k) / Ul{dP e }W[{Pe}] dco e K(u; e )} 



e=i 



(F) 



1 k 
^true l = \ 



p (cos(y>— ip'-ui()) 



(11) 



Notice that so far we did not need to make a choice for the order of the Legendre 
polynomial. One special case of the above equations is in the absence of angular 
disorder and for a fixed degree distribution pk = Sk, c and K{uj) = 5(u). In this case 
we have a c-regular graph in which virtually every spin is living in an identical 
environment and one solution of (11) is ^[{P}] = 5[P((f)) — P*(0)] for some P*(0). 
In the case where the order of the Legendre polynomial is p — 1, so that Li(x) = x 
and the degree distribution is Poisson, our equations reduce identically to those of 
[12] derived via the replica formalism. 

Equation (11) can be solved in a spirit similar to the 'population dynamics' method 
[10,12]. As discussed in the previous section, since Vr[{P}] is a measure over distri- 
butions we can encode each of these distributions using a simple histogram. Once 
a stationary solution for the Vr[{P}] has been obtained one can proceed to evalu- 
ate observables. For example, observables describing magnetization and spin-glass 
order respectively follow from 
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= J { dF}W tIUC [{P}] J d0P(0) cos(A;0) (12) 
gf = J{dF}W tTUC [{P}] [| d0P(0) cos(A;0) 



(13) 



The order parameters (12,13) have the physical meaning 



,.(*) 
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m c' = J™ ^ £ ( cos ( k <Pi)) eq ( 14 ) 

i 

«S*)=kljHS£ (15) 

where (• • -) eg denotes thermal averages (similarly for m^ fc ^ and g^ as averages over 
sin(0)). Generally, the relevant set of observables depends (among other factors) 
on the particular value of p. This is accounted by the dependence of the above on 
the variable k = 1, . . . ,p. One can combine the above observables into a single pair 
e.g. via 

m M = VV^) 2 + (rniy ^) = I( g f + g f) (16) 

For even values of p there is no net magnetization in the system, i.e. rrS^ = 0, 
since the Hamiltonian is then rotationally invariant via a uniform rotation of all 
spins through n. 



3.3 Bifurcation analysis 



To perform a bifurcation analysis, we notice that W^[{P}] = 8(p) P(4>) — sa ^ _ 
isfies (11) for all temperatures. Thus, we can associate this state with the high- 
temperature paramagnetic one (P). If one assumes that bifurcations away from this 
solution occur in a continuous fashion, then we can apply the so-called Guzai ex- 
pansion [12], i.e. consider the perturbation P(<j>) — > ^ + A(0), with / d0A(0) = 
due to normalization. A linear stability analysis then produces the following two 
conditions where, respectively, the first and second moments of W^[{P}] bifurcate: 



l = ^-^max(cos(M)^ (17) 

1 = — max — ^Yvi — (18) 

(k) e>o (jf)y 
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Fig. 1. Phase diagrams of the model (1) for Co = 7r/8 and for different values of p. Left 
figure: Poisson graph with p(k) = e~ c c k /kl. Dashed, solid and dotted lines represent 
the P^F, P— >SG and F^SG transition, respectively. Open circles indicate a change 
in the bifurcating mode. All transitions are second-order. Right figure: c-regular graph; 
markers correspond to integer values of c and lines are guide to the eye. For c = 3 and 
p = 4, 6 the transition is from P^SG. 

where 

/■27T 

= J o # cos(£0) fJ> JL *(<*»*) (19) 

(details on similar calculations can be found in [12]). These two conditions pro- 
vide the critical temperatures where a ferromagnetic (F) and a spin-glass phase 
(SG) appear, respectively. The range of integers over which we maximize (17,18) 
expresses physically the fact that at the moment of bifurcation towards an ordered 
phase the distribution of spin orientations has £ maxima (in fact possible bifur- 
cating modes occur at the Fourier modes cos(£0)). For p = 2 we find that the 
bifurcating mode is given by £ = 2 implying that spins can be ordered in a parallel 
or in an anti-parallel fashion (they are energetically equivalent). For larger values 
of p we find a transition in which the bifurcation mode changes from £ = 2 to 
£ = p. 

Let us now describe the resulting phase diagrams in more detail. For Co = one 
finds a phase diagram with a P^F transition only, while for Co > richer phase 
diagrams occur. In figure 1 we plot phase diagrams for a; = 7r/8 for different values 
of p. The left picture corresponds to Poisson degree distributions p(k) = e~ c c k /k\, 
while in the right picture we show the phase diagram for a regular graph. Larger 
values of p produce sharper P— >F and P— >SG transitions. The critical temperatures 
provided by simulation experiments and the numerical solution of our equations 
(16) are in good agreement with the results of the bifurcation analysis supporting 
the fact that all transitions are second-order (at least for p < 6 where we have 
currently focused). For the F— >SG transition we have assumed that it is given 
by the dotted line. This follows from physical reasoning (absence of re-entrance 
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Fig. 2. Left: The order parameter (16) for the p = 2 model with iD = tt/8 and for 
a Poisson- (dashed) and a regular graph (solid) with (k) = 4. The critical temperatures 
agree well with the results of the bifurcation analysis. Right: The distribution of rod 
orientations P(4>) for p = 2. We observe two peaks at (f> = 0, ir reflecting the energetic 
equivalence of parallel and anti-parallel alignment. Markers in both pictures correspond 
to simulation experiments. 

phenomena [19]). The change in the bifurcation mode is given by the open circles, 
namely for T < T circ i c we find that £ = 2 solves the maximization problem in 
(17,18), whereas for T > T circle we have £ = p. In the left picture of figure 2 we 
show the magnetization order parameter (16) for Poisson and regular random 
graphs of mean connectivity 4. In the Poisson graph transitions are smoother due 
to the variable number of connections per site. Markers correspond to simulation 
experiments of N — 25, 000. In the right picture we show the predicted distribution 
of rod orientations V(4>) = J{dP}W true [{P}]P( y (f)) against simulation experiments 
with V{4>) = 52i{S[(j> — <Pi])eq for a p = 2 model on a Poisson graph with (k) = 5, 
Co = and T = 0.5 and T = 1.5. We observe that P((f>) has two peaks separated by 
7r which reflects the fact that parallel and antiparallel orientations are equivalent. 
Note that we have given an (arbitrary) overall rotation to the spins so that the 
peaks of the distribution are aligned at and ir. Lower temperatures promote order 
and as a result P{4>) is sharper for T/J — 0.5 than for T/J = 1.5. 

Finally, in figure 3 we show the magnetizations m^ 2 \ and for a system 
on a regular lattice with uj — 0. Here we have chosen p = 6 to examine the impact 
of choosing a high-order Legendre polynomial. We observe that at the critical 
temperature T^ t pa 1.16 the magnetisation becomes non-zero; while the other 
two magnetisation order parameters remain zero. This bifurcation occurs at the 
point predicted by (17) which indicates that despite the high p- value the transition 
is continuous. It is interesting to note that the magnetisations and do not 
bifurcate until a lower temperature T^it ~ 1-04 is reached. This transition appears 
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T/J 



Fig. 3. The magnetisations (dashed), (dotted) and (solid), given by equa- 
tion (16), for a system described by a Legendre order p = 6 on a regular lattice and 
with Co = 0. We see that there are two critical temperatures in this system. At the first 
one T^) t w 1.16, which agrees with the prediction of (17), only becomes non-zero. 

For temperature values T < T^ t s=s 1.04 the magnetisations rrS 1 ^ and also become 
non-zero. This latter transition appears to be discontinuous. 

to be discontinuous and, as it is preceded by the transition at T$ t , we have been 
unable to derive it analytically (since in the ordered phase we are required to 
find the distribution P(<f>) numerically). For temperature values between the two 
critical ones the ferromagnetic phase describes a system with 'local' order, i.e. 
relative only to the bifurcating mode (here, the sixth Fourier mode (17)). 



4 Dynamics 

We now turn to the relaxational dynamics of our system. The microscopic dynamics 
of our model is dictated by the Langevin equation 

d , , , dH{4>) 

where r)i(t) represents Gaussian white noise {rji(t}) = and {r)i(t)rjj(t')) = 2T5ijS(t— 
t'). Defining the microscopic state probability p t ((p) = (5(<f> — <f>(t))) (where (• • •) 
represents average over the stochastic process) we can obtain from (20) the corre- 
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sponding Fokker-Planck equation 

^gl^MO (2i) 

We now consider a set of t macroscopic observables which we denote collectively as 
Q(4>) = (fl\((f)), . . . , fle(4>)) and we define accordingly the macroscopic probability 
distribution for these observables: 

v t {n) = J dtptty) 6[Q - n{cf>)} (22) 

From (21) one can derive, through integration by-parts, a Fokker-Planck equation 

for v t (n) 



dt 



Pt{4>) 



N Q 



Pt{4>) 



dH{cj>) 



fc=i 



d 

Ml 



OH dVl k 

d<pi d<pi 



d 2 



h k=1 dQhd^lk 



ft, 



(23) 



which is expressed in terms of the so-called subshell averages 

Sd4>pt(4>)5[n-si(4>)] [•••] 



([•••Da. 



fd<t> Pt (<t>)6[n-n((i>)} 



(24) 



Using (23) we obtain the evolution equations for the macroscopic observables. Let 
us denote = / dQV t (£l) £4- Then integration of (23) gives 



IdH dQ k 

\d(j>i d<f>i 



T- 



(25) 



(assuming that the relevant boundary terms vanish) with the notation ([• • • ])q t — 
J dflp t (ft) ([• • -])$7r Note that (25) while exact (for reasonable choices of fl) is 
not closed as it still depends on the microscopic probability Pt{4>)- To eliminate 
this dependence we now make the equipartitioning assumption which underlies 
dynamical replica theory [16]: we assume that Pt(4>) depends on <fi only through 
f2((/>) so that pt(4>) can be removed completely from the subshell average (i.e. all 
microstates which give the same value of the macroscopic observables are equally 
likely). 



Let us now be more specific and choose as our observables the magnetization and 
energy 
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m 



(k) 



J 



iV 
1 



^ cos(A;0 i ) e = -— L v ( cos (0i - <A? - ^ij)) 



,(*)-_ 
~ N 



in: = — ^sin(A;0j) 



(26) 



(27) 



For notational convenience we will from now suppress the time dependence in 



the above i.e. we will write m 



(k) 

c,t 



m[ k ) and similarly for the other observables. 
Calculating the various derivatives in (25) and inserting into the equations the 
joint spin-field distribution 



D 



(fc) (k) 

1 



(M) 



N 



eedi ax 



5(0-0,) 



mj, ,m> 3 ' ,e,t 



results in the trio of ordinary differential equations 



(2? 



^mW=Tk 2 mW-Jk [ dhd(j)D {k) {k) (h, 0) hsm(k<j)) 
f = Tk 2 m[ k) - Jk j dhd(f)D ( fc) ( fc) (h,<f))hcos(k<f)) 



_d 

~dl 

d 

— m„ 
dt 



d r 
—e = —2Te + J / dhd(f)D ( 
dt J ™c 



(29) 
(30) 
(31) 



The distribution D ( k ) (k) (h, 0) represents the distribution of spins and local 
fields in the microcanonical ensemble for a given value of the magnetization and 
energy. Once we know this at time step t we can predict the values of the observ- 
ables at time step t + At. To calculate this joint distribution our strategy will be to 
associate it to a quantity with the same physical meaning but in another physical 
setting [14]: we consider a canonical ensemble where the temperature and external 
fields will force the system to have the values m^ s t and e t . In this new system the 
Boltzmann distribution is given by 



p{4>) 



l 



-7tf(0-)- 7 £fc>X fc) £.C06(fcfc)+hW £.sin(k&)] 



(fc)' 



Z( % {h { c k \hi k) }) 
while the joint spin-field distribution can be defined as 



(32) 



V(hi, 



d<pdiQ{4>i,4>di)5 



hi 



E sin ( 



dLpjx) 
~ ' dx 



(33) 
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where Q(4>i, <pdi) represents the joint distribution of finding spin % and its neigh- 
bours in the given state. Using the reasoning that follows equation (4) and the 
Bethe approximation (5) we can write for tree-like structures 

Q(<f>h<f>di) ~ e 7 ^ fc > l[ ^ fc)cos(fc0l)+ ^ fe> ^ Yl P {i) {4>t) 

eedi 

(34) 

with p(i\(j)i), as before, representing the cavity distribution of finding spin i in fa 
in the absence of site j and in the symbol ~ we have absorbed the normalization. 
From here, integration with respect to the spin variables of the neighborhood of % 
gives the marginal distribution 

P(<f>i) ~ e 7 ^fc>i [/l ' fc>cos(fc<}!,l)+ ^ fc)sin(fc<}!,l)] f d0 ai e Tj ^eai Lp(cos( ^~^ _a;i£)) ]~J P (i \fa) 

J eedi 

(35) 

and removing a neighbour of % closes this to 



pO')(0.) r~j e 7 Sfc>i[ /l ' fc) cos(fc^)+7^ fc) sin(fc^)] 

x f d<f) dAj e lJ ^^ L P (^-^u)) Y[ pW(^) (36) 

eedi\j 

Observables in this canonical ensemble are given by 

^ = lun^ 1 £ J d0 P(fa) cos(*&) (37) 

i 

^ = Irn^ 1 2 / # P(&) sin(^) (38) 

i 

e=- J5 1 4 51 (L p (cos(fa - fa - Uij)))^ (39) 
^ (»J)e£w 

with the short-hand notation 

/ dfadfaP^{fa)P^\fa) eT J M cos ^-^-^))[. . .] 
^' ~ / d0i P( i )(0 i ) p(i) (0i) e 7ji f( cos ^-^-^')) ^ ^ 

The above equations can be easily re-written in a graph-ensemble form for a given 
degree distribution. 

We see that this is now equivalent to an equivalent static calculation (using a more 
involved Hamiltonian): at each time step t and given the values of observables 
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m c , st and e t , one has to find the values of the temperature 7 and external fields 
h^\h^ that via the definitions (37), (38) and (39) ensure that the equalities 
^(k) = mfj, ^( fe ) = mf], e = et are satisfied. Once the correct values of the 
parameters have been found one can evaluate the joint spin-field distribution (33) 
and this corresponds to D ( k ) 0) (h,(f)). 

While the conceptual formulation of this methodology is relatively straightforward, 
the numerical implementation is a bit more challenging. For each time step, one 
has to solve the inverse problem of finding the correct values of the temperature 7 
and external fields which produce a given magnetization and energy. Thus 

the search space is already (2k + l)-dimensional and computational costs are not 
negligible. The algorithm also requires the convergence of population dynamics and 
we are thus required to work in regions of phase space where only trivial ergodicity 
breaking occurs. For these reasons we have implemented our algorithm in the 
simple scenario of p = 2 (in which case = 0) which allows us to restrict the 
present set of observables to (mf\mf\e). We have considered a 3- regular graph. 
The results are shown in figure 4. We see that the analytic solution compares 
perfectly with simulations for small and very large times, but for intermediate 
times the two clearly deviate. This is due to the assumptions we have made, namely 
equipartitioning in the space of the three observables. We expect that by choosing a 
larger set of macroscopic observables (e.g. incorporating correlation functions) this 
approximation will improve. There is also an error introduced due to the presence 
of loops in the graph; for smaller graphs we would expect a larger deviation due 
to the incorrectness of the Bethe ansatz, while in the truly infinite system these 
will not occur. 



5 Discussion 



In this paper we have studied the Lebwohl-Lasher model, which has a characteristic 
energy landscape of sharp and narrow wells. This model has been studied exten- 
sively in liquid crystal models and our approach here can be seen as the Bethe 
approximation to the finite dimensional problem. We have considered random 
graph lattices. At equilibrium our treatment follows the general methodology de- 
veloped in [11] based on the cavity method within the ergodic (replica-symmetric) 
assumption. We have given bifurcation conditions describing the phase diagram. 
Numerical evaluation of the order parameters shows the transition from ordered to 
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4 1 1 1 1 1 1 1 1 1 1 

U "T) 10 20 30 40 50 
time 

Fig. 4. The evolution of the magnetization mP^ = [m^ + m^]2 for a 3-regular lattice 
and p = 2 from a highly ordered initial state in a heat-bath of temperature T = 0.6 
and Q = 0. The solid line represents the theoretic prediction and markers Langevin 
simulations of N = 40 000 spins. 

paramagnetic phase. This transition, as in Ising spin systems, is sharper for regular 
than for random lattices. Simulation experiments are in excellent agreement with 
the analytic results. 

We have also extended the formalism developed in [13] for the dynamics of spin 
systems on finitely-connected systems into situations where spins are continuous 
variables. Our starting point has been the Langevin equation for the microscopic 
dynamics from which a set of differential equations follows for our chosen set of ob- 
servables. The key to solving this set of evolution equations is the approximations 
of dynamical replica theory [16]. We have contrasted the results of this method- 
ology with simulation experiments. This shows that despite the appealingly clean 
analytic form, there are weaknesses of the formalism that need to be dealt with 
in the future. In particular, the assumption of equipartitioning of the microscopic 
state probability within the subshells is too strong for small sets of observables 
while for larger ones the computational costs start becoming prohibitive. 
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